We calculate perturbatively the tunneling decay rate V 
of the metastable phase in the quantum d = 1 Ising model 
in a skew magnetic field near the coexistence line < h x < 
1, h z — > —0 at T = 0. It is shown that T oscillates in the 
magnetic field h z due to discreteness of the excitation energy 
spectrum. After mapping of the obtained results onto the 
extreme anisotropic d = 2 Ising model at T < T c , we verify 
in the latter model the droplet theory predictions for the free 
energy analytically continued to the metastable phase. We 
find also evidence for the discrete-lattice corrections in this 
metastable phase free energy. 
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It is widely accepted (for review, see 0) that for all 
T < T c , the free energy F(H) of the two-dimensional 
Ising model analytically continued from the positive real 
axis H > into the complex if-plane, has a branch cut 
singularity at the origin. Near the cut drawn along the 
negative real axis H < 0, the imaginary part of the free 
energy is believed to have the following form M 



Im F(e ±J7r \H\)=^B\H\exp {-A/ \H\) 



(1) 



for small \H\. This expression extrapolates to the 
Ising model the results obtained by Langer j| and N.J. 
Giinther et al 0j in the droplet field theory analysis of the 
coarse-grained Ginzburg-Landau model. In the droplet 
theory, the free energy continued to the cut H < is 
interpreted as the free energy of the metastable state 
F ms (H) ee F(e~™ | H |). Langer conjectured ||, that 
Im F ms (H) may be identified (up to a dynamical factor) 
with the metastable phase decay rate provided by the 
thermally activated nucleation. 

The phenomenological droplet theory prediction for 
the amplitude A in @ is @J|,| 



.4 
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(2) 



where M is the spontaneous magnetization, and 
S 2 denotes the square of surface free-energy of the equi- 
librium shape droplet divided by its area. Both E 2 and 
M relate to the equilibrium zero-field state, and are 
known exactly. The quantity E 2 can be calculated by 
use of the Wulff's construction from the exact anisotropic 
surface-tension, as it was shown by Zia and Avron 
and M was obtained by Onsager and Yang 0. The lin- 
ear depending on \H\ prefactor in (|]J) arises in the con- 
tinuum droplet field theory [Q from the contribution of 
the surface excitations (Goldstone modes) of the critical 
droplet. 

Equations ([I]), (|2|) were confirmed by CCA. Giinther, 
Rikvold and Novotny || in numerical constrained trans- 
fer matrix calculations, and by Harris in numerical anal- 
ysis of certain power series [0. However, no analytic mi- 
croscopic evaluation of (|l]) , (g) for the d = 2 Ising model 
did exist. The purpose of thepresent Letter is to perform 
analytic verification of (0) , ^) for the d = 2 Ising model 
in the extreme anisotropic limit. 

We start from the anisotropic Ising model on the 
square lattice. It is characterized by the nearest neighbor 
in-row and in-column coupling constants J\ > 0, J2 > 0, 
the magnetic field H and the inverse temperature (3 = 



l/(fcsT). The extreme anisotropic limit of the model is 
defined as follows B 



r ee cxp(-2/3 J 2 ) -> 0, 



(3) 



h z Ji ->■ 



with constant h x and h z . The transfer matrix of the 
d = 2 Ising model in this limit can be written up to a 
nonsignificant numerical factor as exp(— (3 JiH), where 
H. denotes the quantum spin-1/2 Hamiltonian of the Ising 
chain in a skew magnetic field: 



n = 



N 

E 

n=l 



(4) 



Here a x ' z are the Pauli matrices, N is the number of 
sites in the chain, cyclic boundary conditions are sup- 
posed. The free energy F per the lattice site of the two- 
dimensional Ising model is proportional in limit (|J) to the 
ground state energy E(h x ,h z ) of the quantum Hamilto- 



F = J x lim 



E(h x ,h z ) 



N 



(5) 



Our strategy contains two steps, (i) We calculate per- 
turbatively the tunneling decay rate (that is, the imagi- 
nary part of the energy) of the metastable state in model 
(||) at zero temperature, (ii) By use of mapping @ we 
obtain then the imaginary part of the metastable free en- 
ergy F ms (H) at T < T c for the d = 2 Ising model in limit 

©" 

It should be noted, that models like (|4|) have been 
widely used to describe dynamical properties observed 
in real quasi-one-dimensional Ising ferromagnets [p|,^0[. 
So, the problem outlined in step (i) is by itself of consid- 
erable physical interest. 

In the small field limit /ij « 1, /i 2 « 1 model (Q) was 
studied by Fogedby [jll| . In the free-fermion point h z = 
Hamiltonian (Q) reduces to the form [|l2| 



Ho = H \h, =0 



* d9 



u>{6) ^(0) i/>(0)+ Const, (6) 



where 9 is the quasimomentum, fermionic operators 
if>'(0), ip(6) satisfy the canonical anti-commutational re- 
lations 

{V> f (6>) =2*6(6-0') 
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and 



u(9) = 2 



(l-h x 



4h x sin 2 - 



1/2 



At zero temperature, there is a phase transition point 
h x = 1, which divides ordered (0 < /i x < 1) and disor- 
dered (/la; > 1) phases. Two ferromagnetic ground states 
| 0+) and | 0_) coexist in the interval < h x < 1. They 
are distinguished by the sign of the spontaneous magne- 
tization (0± |< | 0±) = ± M, where M = (1 - hi) 1 / 8 . 

A small negative longitudinal magnetic field h z < 
removes the ground state degeneration. The following 
fermionic representation for the Hamiltonian (Q) is valid 
in the thermodynamic limit N — ► oo 



U = Ho + V + Const, 



(7) 



where Ho denotes the free-fermionic Hamiltonian (^), V 
is given by 



7 = |M^E :ex pf 



(8) 



/7T Jfl _____ 
— exp(^) % /^)(-VW+V t (-^)) 

and ip(9) | + ) = for all 9. We have used the con- 
ventional notation : ... : for the normal ordering with 
respect to the fermionic operators 4 ! (9), ip'(9). Repre- 
sentation (|^), (JsJ) can be obtained from (Q) by applying 
the Jordan- Wigner |D| and duality H transformations. 
In performing the normal ordering of fermionic operators 
in equation we followed to Jimbo et al H|. 

The nonlinear interaction term V in (Q) prevents the 
exact integrability of the model. So, the natural way to 
study model (0) for small | h z |^ is to use a certain 
perturbation expansion. It is clear, however, that the 
straightforward perturbation theory with the zero-order 
Hamiltonian Ho and perturbation V is useless in the con- 
sidered problem. This is due to the fact, that the term V 
contains the long-range interaction Vb between fermions, 
which is given by 



V = V U W ^i=| h z \Mj2- exp(-2 ]T &>,) : (9) 

n^Z j< 7i 

where 

bj = J ^(0)exp(ii0) 

— Vt(0)expHj0) 
Operator Vb is diagonal in the coordinate representation: 



V ° b hn b hn-l- b h I +) 



-^M[7V - 2^(j 2i -hl-x)] C^L-x-^ I °+) ( 10 ) 



i=l 



where j; < jz+i- Since interaction (|1C| ) depends linearly 
on the distance between fermions, it changes the struc- 
ture of the energy spectrum of model (Q) at arbitrary 
small longitudinal magnetic field h z ^ 0. So, to describe 
decay of the metastable vacuum, one should include the 
long range interaction Vb into the zero-order Hamilto- 
nian. This phenomenon is well known in the Stark effect 
p"5| . To describe ionization of an atom by the uniform 
electric field E, one needs to consider the corresponding 
electrostatic energy eEr in a non-perturbative way. 

Accordingly, we subdivide the Hamiltonian (^) into the 
zero-order and interaction parts, as follows 



where 



H = H Q + V 



n =Ho + V + Const, 
V = V - Vb 



(11) 

(12) 
(13) 



The numerical constant in (|T^) in chosen to provide Ho \ 
0+) =0. Since the new zero-order Hamiltonian Ho con- 
serves the number of fermions, its eigenstates can be 
classified by the fermion number. It is clear from dl(j), 
that fermions created by operators b n are just the do- 
main walls dividing the chain into oppositely magnetized 
domains. 

One can easily verify in the small /i^-limit, that the 
metastable vacuum | 0+) decays preferably into a one- 
domain state. We suppose, that this is true also in the 
general case < h x < 1. So, below we shall contract 
the space of considered states to the two-fermion (i.e. 
one-domain) sector. 

Let | 4>i) be the translation invariant two-fermion 
eigenstate of the Hamiltonian Ho ■ In the coordinate rep- 
resentation, the zero-order eigenvalue problem Ho \ <Pi) — 
Ei | (f)) takes the form 

^ 

} J K nn > (j>i{ri) - M \ nh z \ (j>i(n) = —(j>i{n) 

n'eZ 



where 



K i 

- 1 »- r, ii 1 



uj(9) exp[i(rt — n') 



4>i(n) = (0+ | b b n | 4>i), <M~ n ) = -<fii( n ) 

If the energy Ei is small enough, E\ < e, where e <C w(0), 
the wavefunction 4>i(ri) is mainly concentrated far from 
the origin in the classically available region | n \ > oj(0)/(\ 
h z | M). Therefore, we can apply the 'strong coupling 
approximation' JTJj to represent the wavefunction in the 
form 



3 



4>i{n) = ipi(n) - ipi(-n) (14) 
where the function tpi(n) solves the equation 

^ K nn > (fii(n')~ \ h z \ M n (pi(n) = —tpi(n) 

n'eZ 

After the Fourier transform, we obtain 

—<pi(0) exp(m6>) 

where 



ipi(e) = C exp^- 



2\h z \M 
C ={2\h z \ MN)- 1/2 



\f{6) — Ei 0]> (15) 



f(6) = 2 da uj{a) 
Jo 

The 27r-periodicity condition for the function tpi(0) de- 
termines the energy levels Ef. 



El = 1^1 - 2 I h z \ Ml 



(16) 



The normalization constant C in (|T^) is chosen to yield 



•n I W 



where AE = 2 \ h z \ M is the interlevel distance. 

To determine the decay rate T of the metastable vac- 
uum we use the following nonregorous formal procedure. 

(2) 

In the second order correction Ems to the metastable 
vacuum energy we shift the excitation energy levels 
downwards into the complex iJ-plane: 



E^l 



\(0+\V\ ft) 



E, 



(17) 



The width 7 describes phenomenologically the decay rate 
of one-domain states | ^j) . Decay of these states can be 
caused both by term (O) and by other interactions not 
included into the Hamiltonian (Q). As the result, the 
metastable vacuum energy gains the imaginary part 



Im E ms St -7T g(h z ) \(0+\V\ 0,) \\ =0 



(18) 



where 



g(h z ) = Imcot 



/(?r) - Z7T7 

2 \h z I M 



In deriving (|T^) we have extracted from the sum the 



slowly depending on I factor | (0 + | V \ (pi/ 
righthandsidc of ([l7|). 



the 



The metastable vacuum relaxation rate T is deter- 
mined then in the usual way 



T = -2 Im E„ 



(19) 



It is evident from (|T8|), ([l9) that V oscillates in hj 1 with 
the period 2-kM / fjn) . These oscillations become con- 
siderable in the case of coherent tunneling 7 < AE. In 
the limit of noncoherent tunneling 7> AE oscillations 
in h^ 1 vanish and relations ((lj), (|l9] ) transform to the 
Fermi's Gold Rule fWA: 



r = 2tt I (0+ I V I in) I Vo 



(20) 



In the limit h z — > — the matrix element of the interac- 
tion operator can be asymptotically written as 



V\<j>i) 
1 



MN 



2ir W ' d6 



3V2 



(\h z \ MTV) 1/2 



exp 



I f(0o) I 

2 I hr I M 



(21) 



where Oq = i | In h x | is the imaginary zero of the function 
u{9), l>(6 ) = 0. Substitution of ) into @ yields 
finally 

/(go) I 
h, I M 



Im £J ms = --N \h z \M g(h z ) exp \ - 



(22) 

Perhaps, described by ([l9|), ( ]22j ) oscillations in /i z of the 
metastable state decay rate could be observed (indeed, in 
somewhat modified form) in real quasi-one-dimensional 
Ising ferromagnets at very low temperatures. 

Now let us map obtained results to the d 
model. Applying (||) to ( p2| ) we obtain 



Im F ms = B\H I g(H) exp {-A/ \ H |) 



2 Ising 

(23) 



where 



g(H) = Imcot 



Ji [/(*" 



i? I M J 



A 



■h 

M 



f(-ilnh x ) I 



B = —M 
18 



(24) 

(25) 
(26) 



and H — > — 0. These expressions should be compared 
with the droplet theory predictions ([l]), (0). 

First, let us verify, that expressions (|| and ( p5| ) for 
the amplitude A are equivalent. To do this, we need to 
determine the quantity S 2 in limit ([|). 

The droplet equilibrium shajae in the d — 2 Ising model 
is described by the equation 



ai cosh(/3Axi) + a,2 cosh(/3Aa;2) = 1 



(27) 



where x\,x% denote Descartes coordinates of a point on 
the droplet boundary, the scale parameter A determines 
the droplet size, and 



4 



Oi = 



tanh(2/?J 2 



cosh(2/3Ji; 
In the extreme anisotropic limit 

ai S l~2r 2 (l + /i- 2 ), 
and ( |27| ) simplifies to 

Ji 



«2 



_ tanh(2/3J!) 
" cosh(2/3J 2 ) 

a 2 = 4:T 2 /h x 



A 



Integrating in 2:2 this equation we obtain the area of the 
equilibrium shape droplet 5(A) = W/ A 2 , where 

W=-±\f{-i\nh x )\ 

P 

It follows from the Wulff's theorem ||, that the surface 
energy S(A) also can be expressed in W: S(A) = 2VF/A. 
Therefore, E 2 = AW, and 

A = — = — I f(-ilnh x ) I 
8M M 71 



in exact agreement with ( 2E 

Further, expression (^3) differs from (jl]) by the oscil- 
lating factor g{H). We interpret this factor as the correc- 
tion coursed by the discrete-lattice effects. Those may be 
significant at low temperatures in the presence of strong 
anisotropy p"§| ] . The following observation supports such 
interpretation. 

At low temperatures (h x — > 0), the factor g(H) can be 
written as 



g(H) = Imcot <^ tt • 2x 1 (H) 



iir 7 Ji 
2 \H\M 



(28) 



where 2xi(H) — 2Jx/(M\H\) is the continuum nucle- 
ation theory value of the critical droplet diameter in the 
xi-direction. Maximum points in ( p8| ) just correspond to 
discrete values of this diameter. 

The oscillatory factor g(H) contains parameter 7 
which remains undetermined in the present incomplete 
theory. One would expect, however, that in the critical 
region (h x — > I) parameter 7 is large enough, so that 
g(H) = 1, and oscillations in ( p3| ) vanish. Really, in 
this limit spectrum (|l6| ) becomes continuous, and T can 
be obtained directly from ( pp| ) without referring to (17), 

In the critical region expression ( p6[) agrees with our 
previous result |l^] . 
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